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INTRODUCTION 

The  mechanical  response  and  the  failure  of  composite  solid 
propellants  are  known  to  be  related  to  the  formation  and  growth  of 
vacuoles  on  the  microscopic  or  macroscopic  scale.  Failures  in  com- 
posite materials  such  as  coiiq>oslte  propellants  originate  at  the  filler 
particle  or  binder  molecular  level.  These  microscopic  failures, 
whether  adhesive  failures  between  the  binder  and  filler  particles  or 
cohesive  failures  in  the  binder,  result  in  vacuole  formation.  The 
cumulative  effect  on  all  vacuole  formation  and  growth  is  observed  on 
the  macroscopic  scale  as  strain  dilatation.  The  strain  dilatation  in 
composite  propellants  varies  with  binder  Lype,  binder  formulation, 
filler  particles,  and  presumably  filler  particle  size  distribution. 

The  strain  dilatation  caused  by  vacuole  growth  and  formation  results 
in  a nonlinear  stress-strain  behavior  for  these  materials. 


The  experimental  investigations  of  Leemlng  et  al.  (1), 
Bornsteln  (2),  and  Martin  (3)  emphasize  the  fact  that  the  state-of-the- 
art  linear  viscoelastic  theory  does  not  predict  propellant  grain 
stresses  and  strains  under  realistic  loading  conditions.  The  most 
significant  portion  of  the  error  is  due  to  the  constitutive  nonlinear- 
ity of  the  material.  Linear  thermorheological  simple  viscoelastic 
theory  has  certain  important  limitations;  however,  it  is  still  the  only 
constitutive  theory  currently  used  in  most  analyses  because  no  other 
constitutive  equations  previously  existed  that  are  as  easily  incorpo- 
rated in  the  analyses. 
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There  has  been  considerable  progress  in  the  development  of 
nonlinear  constitutive  theories  and  methods  of  approximating  specific 
nonlinear  viscoelastic  behavior  in  the  last  few  years.  There  will  be 
no  attempt  in  this  paper  to  review  the  extensive  literature  on  the 
general  theories  of  nonlinear  viscoelasticity. 

BACKGROUND  AND  THEORY  . 

Farris  (4)  demonstrated  that  the  mechanical  response  and  the 
failure  in  solid  propellants  are  related  to  the  formation  and  growth 
of  vacuoles  within  the  composite  material  and  that  this  phenomenon  is 
directly  related  to  the  macrcscopically  observed  strain  dilatation. 

The  mechanical  properties,  particularly  Poisson's  ratio,  are  strongly 
dependent  on  the  dilatatlonal  behavior  cf  these  materials.  Strain 
dilatation  measurements,  even  for  the  uniaxial  tensile  tests,  are 
cumbersome,  tedious,  and  time  ccntumlng.  Therefore,  a relationship 
that  would  permit  the  determination  of  strain  dilatation  and  Poisson's 
ratio  as  a function  of  strain  utilizing  only  the  stress -strain  data 
would  be  most  desirable.  Insight  into  the  factors  which  govern  the 
nonlinear  response  and  failure  of  composite  propellants  can  be  obtained 
utilizing  a model  for  dewetting,  vacuole  formation,  and  growth  of 
vacuoles  which  is  based  on  assumptions  regarding  the  micros true tural 
behavior  of  the  system.  The  model  can  then  be  used  to  develop  an 
expression  for  dilatation  as  a function  of  stress  and  strain  and  to 
conqjare  the  predicted  dilatation  with  experimental  measurements  of 
dilatation.  The  degree  of  agreement  between  the  predicted  and  experi- 
mental values  of  dilatation  gives  an  indication  of  the  validity  of  the 
model.  This  approach  was  used  by  Farris  (5,6),  Fedors  and  Landel  (7), 
and  others  (8,9). 

The  present  state  of  the  art  of  predicting  dilatation  as  a 
function  of  strain  is  limited  because  there  is  very  littlo  quantita- 
tive Information  available  on  the  experimentally  determined  dilatation 
under  conditions  similar  to  those  encountered  with  solid  propellant 
rocket  motors.  Various  models  describing  this  phenomenon  have 
appeared  in  the  literature.  In  developing  such  a model,  attention  is 
focused  on  a rigid  filler  particle  contained  in  an  elastomeric  matrix. 
As  the  material  is  strained  above  some  critical  value,  a vacuole  is 
formed  about  the  filler  particle  due  to  the  internal  failure  of  the 
composite.  The  source  of  vacuole  formation  may  be  either  adhesive 
failure  of  the  filler-binder  bond  or  cohesive  failure  in  the  binder 
near  the  filler  particle  (10,11).  As  the  strain  in  the  material 
Increases,  the  vacuoles  continue  to  Increase  in  size  and  number.  The 
shape  and  instantaneous  behavior  of  these  vacuoles  seem  to  follow 
these  assumptions: 
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1)  Vacuoles  from  arbitrarily  at  any  magnitude  of  strain 
above  the  critical  value. 

2)  Each  vacuole  behaves  as  an  ellipsoid  of  revolution  with 
the  minor  axis  determined  by  the  diameter  of  the 
enclosed  filler  particle. 

3)  The  major  axis  of  the  ellipsoidal  vacuole  increases 

• linearly  witlk  strain  at  a rate  proportional  to  the  size 

of  the  filler  particle  it  contains. 

Various  constitutive  theories  exist  and  experimental  investi- 
gations can  be  designed  to  produce  specific  stress-strain  equations. 

In  this  section  a constitutive  relationship  for  uniaxial  tensile  data 
is  discussed  and  the  similarity  with  Farris'  equation  is  noted  (4). 

Martin  (12)  used  the  constitutive  relationship 

a(e,T,t)  - {F(e,T,t)e‘^"^^®^}6  . (1) 

Equation  (1)  has  been  verified  experimentally. 

Starting  with  Farris'  equation  and  defining  F ■ o/s  as  the 
secant  modulus,  as  was  used  by  Martin,  the  following  relationship  is 
obtained : 


F(e) 


(2) 


\diere  is  the  initial  modulus  and  FCe)  is  meant  to  indicate  F at  a 

given  strain.  This  value  will  be  indicated  simply  by  F in  the 
remainder  of  this  work.  After  rearrangement  one  obtains: 


JL  , -3(AV/V)l/e 


(3) 


Equation  (3)  indicates  that  the  nonlinearity  of  the  stress-strain  curve 
(of  i«hich  the  vacuole  formation  and  growth  is  considered  to  be  the 
major  contributor)  shoitld  be  described  by  the  ratio  of  F/E^.  Taking 


the  natural  logarithm  of  both  sides  of  Equation  (3)  results  in  Equa- 
tion (4): 


<4) 


5.  -WWi.c' 
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Since  3 is  assumed  to  be  Independent  of  strain  magnitude,  Equation  (3) 
suggests  that  a plot  of  {ein(E^/F))  versus  strain  should  have  a simi- 
lar shape  to  the  dilatation-strain  curve.  Figure  1 shoves  plots  of 
{6in(Ej^/F))  and  (AV/V)  versus  strain.  The  similarity  in  the  shapes  of 


Figure  1.  Total  Nonlinearity  and  Dilatation  Versus 
Strain  for  a Granular  Filled  Elastomer. 


the  two  curves  indicates  that  the  total  nonlinearity,  {€in(Ej^/F)  ),  in 

the  materials'  behavior  may  be  calculated  from  an  equation  similar  to 
that  used  by  Farris  (4).  The  same  reasoning  used  by  Farris  results  in 
the  expression 


f(€) 


e-n2/2 

dndn  + y e 


(5) 


where 


f(€)  * total  nonlinearity  of  the  materials'  behavior 


S 


max 
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7^  = initial  slope  of  the  {ein(Ej^/F)  ) versus  strain  curve 
7 ■ final  slope  of  the  {ein(E,/F))  versus  strain  curve 

lUflX  A 

1 " the  strain  magnitude  at  the  intersection  of  the  initial 
asymptote  of  slope  7 and  the  final  asymptote  of  slope 

7 ° 

. 'max 

A « the  magnitude  of  (€in(Ej/F))  at 

B B the  magnitude  of  the  ordinate  at  the  Intersection  of  the 
two  asymptotes 

n » (e  - €)/S. 

Equation  (5>  then  would  permit  one  to  estimate  the  total  nonlinearity 
cf  the  propellant  once  the  parameters  have  been  determined.  It  would 
be  more  useful,  however,  to  obtain  a quantity  representing  the  pro- 
pellant's nonlinear  behavior  that  could  be  readily  utilized  in  the 
current  linear  viscoelastic  programs  to  account  for  the  nonlinear 
response  of  the  material.  With  this  In  mind,  consider  Equation  (3), 
which  after  rearranging  becomes 


e{f)  - 


(6) 


Recalling  that  ^ and  E^  are  independent  of  strain  and  differentiating 
both  sides  with  respect  to  strain,  one  has 


(7) 


where  E the  Instantaneous  slope  of  the  stress-strain  curve.  Equa- 
tion (7)  now  contains  two  unknowns,  ^ and  d/de(AV/V  ).  Now,  consider 
I o 

the  change  of  dilatation  with  respect  to  strain.  According  to  Farris' 

data  (5),  the  quantity  d/de(AV/V  ) obtains  a maximum  somewhere  between 

o 

the  points  of  maximum  stress  and  rupture  stress  and  remains  constant 
until  failure  occurs.  Farris  proposed  the  relationship 


- c V 


» 


(8) 
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where  ” volume  fraction  of  dewetted  solids*  Experimental  data 

were  presented  for  systems  that  tend  to  dewet  completely.  These  data 
gave  excellent  agreement  with  the  model  represented  by  Equation  (8). 
The  maximum  slope  of  the  dilatation-* strain  curve  is  shown  to  be 
directly  proportional  to  the  volume  fraction  of  filler  particles  in 
the  material.  For  materials  that  tend  to  dewet  completely,  the  pro- 
portionality factor  is  unity  (c  ■ 1) ; therefore, 


(9) 


where  is  the  total  volume  fraction  of  filler  particles  in  the 
system. 


With  these  observations  noted,  the  following  simplifying 
assumptions  will  be  made  that  will  enable  the  determination  of  the 
factor  P without  having  to  conduct  the  dilatation-strain  measurements: 

1)  The  maximum  slope  of  the  dilatation- strain  curve  is 
obtained  at  a strain  magnitude  equal  to  halfway  between 
the  point  of  maximum  stress  and  breaking  stress. 

2)  All  composite  propellants  are  totally  dewetted  in  the 
region  of  the  sample  where  falluie  occurs. 

The  assumption  that  all  filler  particles  in  composite  pro- 
pellants are  totally  dewetted  in  the  region  of  failure  is  not  too 
unrealistic  for  the  following  reasons.  In  composite  propellants  vdiere 
the  bond  between  the  filler  particles  and  the  binder  is  stronger  than 
the  cohesive  strength  of  the  binder,  vacuoles  Initiate  in  the  binder 
and  propagate  to  the  surface  of  the  filler  particle  prior  to  failure. 
In  these  materials  there  is  probably  a large  strain  gradient  in  the 
region  of  failure  and  there  would  be  a high  probability  of  totally 
dewetted  filler  particles  In  this  region.  One  reason  that  Equation  (9) 
does  not  appear  to  be  true  for  all  propellants  is  that  the  measured 
dilatation,  whether  localized  or  not, Is  averaged  over  the  total  volume 
of  the  sample.  If  the  localized  band  of  dewetted  filler  particles  Is 
narrow,  the  measured  dilatation  would  appear  to  be  ouch  smaller  than 
the  actual  dilatation  In  the  region  of  failure.  The  maximum  slope  of 
the  dllatatlon-straln  curve  Is  related  to  the  dewetted  filler  parti- 
cles, but  the  dewetted  filler  particles  may  be  contained  In  a narrow 
band  around  the  failure  region. 
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The  assumption  chat  the  maximum  slope  of  the  dilatation- 

strain  curve  is  obtained  at  a strain  magnitude  halfway  between  the 

points  of  maximim  stress  and  breaking  stress  is  approximately  correct 

according  to  Farris'  data.  The  purpose  of  this  assumption  is  to  fix 

the  point  on  the  stress-strain  curve  for  the  calculation  of  the 

coefficient  0.  One  reason  that  failure  points  are  not  recoonended  for 

this  calculation  is  the  variation  in  ways  of  selecting  failure  points. 

The  use  of  the  stress  ana  strain  values  at  £ •■  (e  +£.1/2  should 

m D 

minimize  the  variations  in  methods  of  selecting  failure  points  as  a 
source  of  error  in  the  calculation  of  3. 


Substituting  Equation  (9)  into  Equation  (7)  yields  the 
following  results  after  rearranging: 


P 


(10) 


to  be  evaluated  at  € ■ (£  + e.)/2. 
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With  ^ known,  one  can  then  utilize  Equation  (6)  to  determine 
the  dilatation  as  a function  of  stress  and  strain.  The  resulting 
expression  is  as  follows: 


This  equation  was  utilized  to  calculate  the  apparent  volume 
change  for  two  composite  materials  where  the  experimental  dilatation- 
strain  data  were  available.  The  calculated  and  experimentally  measured 
dilatation- strain  data  were  in  excellent  agreement. 

To  more  readily  utilize  the  volume  change  data  obtained  in 
this  way  for  a nonlinear  viscoelastic  stress  analysis,  the  following 
method  is  proposed.  It  would  be  more  desirable  to  reflect  the  pro- 
pellant nonlinearity  in  an  equivalent  Poisson's  ratio,  which  would  be 
allowed  to  vary  with  strain  magnitude,  strain  rate,  and  temperature. 

The  variable  Poisson's  ratio  could  then  be  Incorporated  into  the  cur- 
rent linear  viscoelastic  stress  analysis  program  to  reflect  the  pro- 
pellant's nonlinear  behavior. 

To  derive  an  expression  of  Poisson's  ratio  as  a function  of 
stress  and  strain,  the  volume  change  is  approximated  by  the  first 
strain  Invariant  neglecting  higher  order  terms  involving  strain. 
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vhere  62*  and  are  the  strains  in  the  principal  directions  and 
are  related  as  follows  in  the  uniaxial  condition: 


-ve. 


-ve 


(13) 


where  v is  Poisson's  ratio.  Substituting  these  relationships  into 
Equation  (11)  results  in  the  following 


(14) 


Substituting  Equation  (14)  into  Equation  (6)  and  rearranging,  one 
obtains  the  relationship: 


V 


1 


(15) 


Also,  substituting  Equation  (14)  into  Farris'  equation  yields  the 
following  constitutive  relationship  for  the  uniaxial  nonlinear 
behavior: 


0 ■ Ej^ee 


-P(l-2v) 


(16) 


For  the  nonlinear  viscoelastic  treatment,  the  initial  slope  modulus, 
£^,  is  allowed  to  vary  with  strain  rate  and  temperature;  the  non- 
linearity coefficient,  p,  is  allowed  to  vary  with  strain  rate  and 
temperature;  and  the  Poisson's  ratio,  v,  is  allowed  to  vary  with 
strain  rate,  temperature,  and  strain  magnitude. 


EXPERIMENTAL  VERIFICATION 


The  model  used  in  deriving  the  expressions  for  dilatation 
and  Poisson's  ratio  as  a function  of  stress  and  strain  was  based  on 
certain  assumptions  regarding  the  micros true tural  behavior  of  com- 
posite materials  such  as  composite  propellants.  The  validity  of  uhe 
assumptions  and  the  model  is  indicated  by  good  agreement  of  calculated 
and  experimental  stress  and  dilatation  data. 
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Figures  2,  3,  snd  4 are  plots  of  stress  and  dilatation 
versus  strain  for  a granular  filled  elastomer.  The  experimental  data 
were  obtained  from  Farris 'work  (5).  The  experimental  data  in  these 
figures  were  obtained  at  ZS^C  (75*F)  and  at  strain  rates  of  6.66, 

66.6,  and  666X/min.  The  plotted  points  represent  the  data  read  from 
Figure  6 in  Farris'  paper.  The  solid  lines  indicate  the  values  cal- 
culated using  Equation  (16)  for  stress  and  Equation  (11)  for  dilata- 
tion! Figures  2,  3,  and  4 indicate  good  agreement  between  calculated 
and  experimental  values  of  stress  and  dilatation  at  all  strain  rates. 
The  dilatation-strain  data  appear  to  be  shifted  by  a constant  factor 
along  the  strain  axis  in  some  cases.  This  would  Indicate  a slight 
error  in  the  dilatation  measurements,  which  could  be  attribuned  to 
localized  dewetting  or  vacuole  formation  as  discussed  previously. 
Equation  (16)  Is  shown  to  yield  excellent  agreement  with  experimental 
stress  values  at  all  strain  rates.  In  Figure  4,  Equation  (11)  is 
sho%m  to  also  give  excellent  agreement  with  the  dilatation  experimen- 
tal data  at  a 8tr;>in  rate  of  666%/min  up  to  a strain  magnitude  of  35%. 
At  strain  rftagnitudes  in  excess  of  35%,  the  experimental  dilatation 
data  fall  below  the  predicted  values.  VRiile  vacuoles  may  form  arbi- 
trarily at  all  strain  magnitudes  above  a critical  value,  it  is  reason- 
able to  assume  that  the  weakest  point  in  the  composite  toaterial  will 
be  more  susceptible  to  vacuole  formations  and  should  contain  a larger 
population  of  vacuoles  than  the  remainder  of  the  material.  The  forma- 
tion of  vacuoles  to  relieve  the  large  strain  gradient  that  would  be 
present  at  the  weakest  point  in  the  composite  requires  a finite  time 
in  which  to  react.  If  the  composite  material  is  strained  at  a rate 
faster  than  the  material  can  respond  through  mlcrostructural  failure 
and  vacuole  formation,  the  dilatation  would  appear  to  be  uniform 
throughout  the  material  until  conditions  were  such  that  the  rate  of 
dilatation  at  the  weakest  point  in  the  material  became  large  enough  to 
promote  the  formation  of  a localized  band  of  a large  population  of 
vacuoles.  The  data  presented  in  Figure  4 indicate  that  at  a strain 
rate  of  66C%/mln, near  uniform  dilatation  is  present  up  to  a strain 
magnitude  of  approximately  35%.  At  this  strain  a localized  band  of  a 
large  population  of  vacuoles  developed  at  the  weakest  point  in  the 
composlce  material,  which  then  began  to  have  a strong  influence  on  the 
material's  behavior.  As  the  material  is  strained  in  excess  of  35%, 
the  volume  increase  in  the  vacuoles  within  the  established  localized 
failure  region  is  more  than  the  measured  average  dilatation.  This 
would  explain  why  the  experimental  dilatation  data  presented  in 
Figure  4 became  less  than  the  calculated  values  as  the  strain  magni- 
tude Increased  to  exceed  35%.  At  the  lower  strain  rites  of  6.66  and 
66.6%/min,  the  localized  failure  region  was  established  in  the  initial 
portion  of  the  test  (e  < 0.10),  thus  explaining  why  experimental 
dilatation  deviated  slightly  from  the  calculated  dilatation  for  the 
duration  of  the  test. 
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Figure  2.  Dllatetlon  end  Stress  Versus  Strain 
for  a Granular  Filled  Elastomer, 


Figure  3.  Dilatation  and  Stress  Versus  Strain 
for  a Granular  Filled  Elastomer. 


DILATATION  (AV/V)  DILATATION  (AV/Vt 
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Figure  4.  DllaCetion  and  Stress  Versus  Strain 
for  a Granular  Filled  Elastomer. 


Figures  2,  5,  and  6 are  plots  of  stress  and  dilatation 
versus  strain.  The  experimental  data  In  these  figures  were  obtained 
from  Figure  3 In  Farris'  paper  (5).  These  data  Indicate  the  stress 
and  dilatation  behavior  of  a typical  granular  filled  elastomer  at  a 
strain  rate  of  6.66%/mln  at  2 5" , 4.4*,  -18"C  (77®,  40®,  and  0®F) . 
The  data  shown  or.  Figures  2,  5,  and  6 Indicate  excellent  agreement 
between  values  of  stress  calculated  by  Equation  (15)  and  the  experi- 
mental values  at  all  temperatures.  The  behavior  of  the  dilatation 
strain  data  Is  similar  to  that  shown  on  Figures  2,  3,  and  4 with  the 
experimental  values  deviating  slightly  from  the  values  calculated  by 
Equation  (11).  The  same  reasoning  as  used  previously  would  also 
explain  the  deviation  of  the  experimental  dilatation  from  that  cal- 
culated by  Equation  (H) , as  Indicated  on  Figures  2,  5,  and  6. 

Normally,  nominal  stress  values  were  corrected  to  true 
stress  by  multiplying  the  nominal  stress  value  by  the  extension 
ratio,  This  correction  was  originally  derived  based  on  the 
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Figure  5.  Dilatation  and  Stress  Versus  Strain 
for  a Granular  Filled  Elastomer. 
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Figure  6.  Dilatation  and  Stress  Versus  Strain 
for  a Granular  Filled  Elastome/. 
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assumption  that  there  Is  no  volume  change  within  the  sample  for  the 
duration  of  the  test.  Ihe  relationship  considering  volume  change  is 
derived  as  follows: 


1 


(17) 


where  ■ 1 + 


Therefore,  rearranging  Equation  (17),  one  obtains 


(18) 


Substituting  this  expression  into  the  equation  o * o A /A,  one 
obtains  the  following: 


n 


c 


a 

o 1 


1 + 


AV 


(19) 


Therefore,  the  true  stress  or  is  lower  than  would  be  obtained  by  the 
normal  correction  procedures.  Equation  (19)  reduces  to  the  equation 
normally  used  as  AV  approaches  zero. 


CONCLUblONS 

A nonlinear  constitutive  relationship  is  developed  for 
uniaxial  stress-strain  conditions  that  are  shown  to  fit  experimental 
data  over  a wide  range  of  strain  rates  and  temperatures.  The  total 
stress-strain  nonlinearity  of  particulate  filled  composites  may  be 
represented  by  statistical  parameters  similar  to  those  used  by  Farris 
to  describe -the  dilatation-strain  behavior By  using  Equations  (10) 
and  (11)  ,'%he  dilatatlonal  bSEiavibr  of  composite  propellants  may  be 
predicted  fifom  the  stress-strain  data  and  the  total  filler  content. 

The  relationships  derived  in  this  report  will  enable  one  to  investi- 
gate the  nonlinear  behavior  at  strain-rates  and  temperatures ■ vdiere  it 
is  very  difficult,  if  not  impossible,  to  obtain  measured  dilatation 
data. ^ 

\ 

The  relatio''anip  represented  by  Equation  (15)  may  be  uti- 
lized, along  with  experimental  stress-strain  data  obtained  at  differ- 
ent strain- rates  and  temperatures,  to  determine  the  dependency  of  an 
equivalent  Poisson's  ratio  on  the  strain-rate,  temperature,  and  strain 
magnitude.  The  Poisson's  ratio  thus  determined  may  then  be  utilized 
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in  currently  available  computer  programs  to  describe  transient  thermo- 
vlscoelastlc  behavior  with  the  consideration  of  the  materials’  non- 
linear response. 

The  relation  normally  used  for  obtaining  corrected  stress 
data  by  multiplying  the  nominal  values  by  the  extension  ratio  would 
yield  values  of  stress  and  strain  at  the  maximum  corrected  stress, 
\^ich  were  significantly  larger  than  the  true  stress-strain  values 
would  Indicate.  These  errors  could  conceivably  result  In  the  wrong 
conclusions  as  to  the  structural  integrity  of  marginal  designs;  I.A., 
one  might  predict  a safe  design  utilizing  stress-strain  allowables 
based  on  the  normal  correction  method  when,  In  reality,  the  true 
stress-strain  allowables  would  indicate  the  design  would  be  expected 
to  fall.  It  is  therefore  recommended  that  stress-strain  allowables 
for  use  In  future  failure  analysis  be  based  cn  true  stress  as  defined 
by  Equation  (18). 
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